ON THE COEFFICIENTS OF 
THE ASYMPTOTIC EXPANSION OF n\ 
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Abstract. Applying a theorem of Howard for a formula recently proved by Brassesco and 
Mendez, we derive new simple explicit formulas for the coefficients of the asymptotic expan- 
sion to the sequence of factorials. To our knowledge no explicit formula containing only the four 
basic operations was known until now. 



1. Introduction 

It is well known that the factorial of a positive integer n has the asymptotic expansion 



(1.1) n!~n n e-™v^53^, 



n 

k>0 



known as the Stirling's formula (see, e.g., [IIEIII]). The coefficients a k in this series are usually 
called the Stirling coefficients [TJ |6] and can be computed from the sequence b k defined by the 
recurrence relation 

(1.2) b k = -j-Lj L k _x -J^jbjbk-j+^j ,60 = 61 = 1, 

as a k = (2k + l)!!62fc+i 13 II]- It was pointed out by Paris and Kaminski [6] that "There is no known 
closed-form representation for the Stirling coefficients" . However there is a closed-form expression 
that involves combinatorial quantities due to Comtet [5]: 

2k 

(1.3) ak =j2(-iy 

3=0 

where (I3 (p, q) is the number of permutations of p with q permutation cycles all of which are > 3. 
Brassesco and Mendez proved in a recent paper [7] that 

(1 a\ 1 , j S 3 (2k + 2j,j) 

afc = g (_1) 2^ { k +j) r 

where S3 (p, q) denotes the 3-associated Stirling numbers of the second kind. We show that the 
Stirling coefficients a k can be expressed in terms of the conventional Stirling numbers of the second 
kind. A corollary of this result is that the Stirling coefficients have a representation that involves 
only the four basic operations, i.e., an explicit, exact expression. 



y r n j d 3 (2k + 2j,j) 
2^ > 2 k +i{k + j)\ 
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2. The formulas for coefficients 
One of our main results is the following: 
Theorem 2.1. The Stirling coefficients have a representation of the form 

_(2k)\™(k + i-l/2\(3k+l/2\*(i\ l S(2k + l + J ,j) 

where S (p, q) denotes the Stirling numbers of the second kind. 
From the explicit formula 



\2k+i+j 



q ' 1=0 v / 

we immediately obtain our second main result. 

Corollary 2.1. The Stirling coefficients have an exact representation of the form 

Note that this formula contains only the four basic operations. To prove Theorem 12.11 we need 
some concepts. Let r > and a r ^ 0, let F (x) = Y^j> r a j x ''/^- ^ e a f° rma l power series. The 
potential polynomials F„ in the variable z are defined by the exponential generating function 

v v ' ' n>0 

For r > 1, the exponential Bell polynomials B n ^ (0, . . . , 0, a r , a r +i> • ■ •) in a n infinite number of 
variables a r , a r +i, . . . can be defined by 

Ex n 
B ni< (0 ) ...,0,o r ,o r+ i,.-.)-r- 

n>0 

The following theorem is due to Howard [2]. 

Theorem 2.2. If Fn * s defined by (|2.3[) and B n i is defined by (j2.4[) , then 
c_i , /z + i — 1 \ ( z + n\ / r! \ 1 

(2.5) ^=E(-ir( , ^^(o,..^,^^,..). 

Now we prove Theorem 12. 11 
Proof of Theorem \2.1\ Brassesco and Mendez showed that if 



U + 2)! 



then 

(2.7) a , = _L^(f; ^)«». 
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where d k / denotes the fcth derivative of a function /. Define the polynomials Gn^ in the variable 
z by the following exponential generating function: 

(2.8) V-vo"* 



2 e x — x — 1 / 4-^ 3 ? ! 

j>0 



Inserting z = J into this expression gives 



On the other hand we have by series expansion 

(2.10) G-^( X )=J2V(G-^)(0) X -L. 

Equating the coefficients in (|2.9I) and (12.10[) gives 

cP(G- 2J ^) (0) = gS-L G H). 
Now by comparing this with (j2.7[) yields 
(2-11) « fc = ^!^ +i \ 

Putting r = 2 an a r = a r+ \ = . . . = 1 into the formal power series i* 1 (a;) = X)j> r a i a;J '/ gi yes 
i* 1 (x) = — x — 1. And therefore the generated potential polynomials are 



2k + l 



r Qi z )_ 

e x — x — 1 J \2 e x — x — 1 / ^— ( 3 j ! 



j>0 



According to Howard's theorem we find 



<w <*> = 1 <-»* e + r l ) (r>^*~« <«. u, . . .>. 
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Now we derive an expression for the exponential Bell polynomials B n ^ (0, 1,1,.. .) in terms of the 
Stirling numbers of the second kind: 

i\ Bn ti (0, 1, 1, . . .) — r = (F (x)Y = (e* -x- If 

n>0 

=£(;)(-!>'-' ^(ett 

1>1 ' I j=0 V/ \l>l 



j=0 n>0 



l 

X" 



EEC) (-'r J ^ ( n,i) , 

n>0j=0 V/ 



n\ -A / A i_j S{n-i+ j,j) I x r ' 



«E ?E <-ir^ 



. i\ \j J (n — i + ?)! I n\ 



Hence 
(2.13) 

Thus we obtain 



(2,4) «'=e ! r_* 2 *« ! E ' hi 



Substituting z = k + 1/2 and n = 2k into this expression yields 



2 A- 



(2.i5, G <r » - e r + V 1/2 ) f » + : '? > : ^ (2fc+i+ "> 



„ 2fe-i y~ v "' v '^Vi7 v ~ y J ' (2k + i+j)\ 



hence by (|2.11[) we finally have 



<-) « - §§f | C +i - 1/2 ) C£!f)*g G) '-""^^f 1 - 

This completes the proof of the theorem. □ 
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